Effective temperature in driven vortex lattices with random pinning 
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We study numerically correlation and response functions in non-equilibrium driven vortex lattices 
with random pinning. From a generalized fluctuation-dissipation relation we calculate an effective 
transverse temperature in the fluid moving phase. We find that the effective temperature decreases 
with increasing driving force and becomes equal to the equilibrium melting temperature when the 
dynamic transverse freezing occurs. We also discuss how the effective temperature can be measured 
experimentally from a generalized Kubo formula. 
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Whether and how can one extend thermodynamic con- 
cepts to nonequilibrium systems is a very important chal- 
lenge in theoretical physics. Many definitions of nonequi- 
librium temperatures have been proposed in different 
contexts but it has been rarely checked if they conform 
with the expected properties of a temperature. 

Cugliandolo, Kurchan and Peliti Q have introduced 
the notion of time-scale dependent "effective temper- 
atures" T a ff from a modification of the fluctuation- 
dissipation theorem (FDT) in slowly evolving out of equi- 
librium systems. T eil is defined from the slope of the 
parametric plot of the integrated response against the 
correlation function of a given pair of observables when 
the latter is bounded or, equivalently, it is the inverse of 
twice the slope of the parametric plot of the integrated 
response against the displacement when the correlation 
is unbounded. This definition yields a bona fide tempera- 
ture in the thermodynamic sense since it can be measured 
with a thermometer, it controls the direction of heat flow 
for a given time scale and it satisfies a zero-th law within 
each time scale. T eff was found analytically in mean- 
field glassy models [|| and it was successfully studied in 
structural and spin glasses, both numerically and ex- 
perimentally p|, in granular matter and in weakly 
driven sheared fluids R || . 

In their study of driven vortex lattices in type II su- 
perconductors, Koshelev and Vinokur || have defined 
a "shaking" temperature T sll from the fluctuating force 
felt by a vortex configuration moving in a random pin- 
ning potential landscape. This lead to the prediction of 
a dynamic phase transition between a liquid-like phase 
of vortices moving at low driving forces and a crystalline 
vortex lattice moving at large forces, when T sh equals 
the equilibrium melting temperature of the vortex system 
H [n|. However, later work pd| O , Q has shown that 
the perturbation theory used in p] breaks down and that 
the vortex phase at large velocities can be an anisotropic 
transverse glass instead of a crystal. In spite of this, the 



shaking temperature introduced in |9(] has been a use- 
ful qualitative concept, at least phenomenologically. In- 
deed, the dynamic transitions and moving vortex phases 
discussed inj 
imentally Il4 



. [To|. [l2, 13 have been observed exper- 
and in numerical simulations 15, fL6|. 



In this Letter we apply the definition of T eff based 
on the modifications of the FDT to driven vortex lattices 
above the critical force and within the fluid moving phase. 
We compare our results with the shaking temperature of 
H and discuss how to obtain T eli experimentally from 
measurements of transverse voltage noise and transverse 
resistance p7| . 

The equation of motion of a vortex in position R; is: 
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where Rij = 

Rip — 



IRi 



Rj| is the distance between vortices 
^,J, Kip = |Rj — R p | is the distance between the vor- 
tex i and a pinning site at R p , rj is the Bardeen-Stephen 
friction, and F = ^£0 J x z is the driving force due to a 
uniform current density J. The effect of a thermal bath 
at temperature T is given by the stochastic force Ci(^)j 
satisfying (£<*(*)) = and (£<* (t)CjV)> = 2 V k B T8(t - 
t^SijSfj,^' , where (. . .) denotes average over the ensemble 
of d . We model a 2D thin film superconductor of thick- 
ness d and size L by considering a logarithmic vortex- 
vortex interaction potential: U v (r) = — A v ln(r/A), with 
A v = $g/8?rA and A = 2X 2 /d > L (l6). The vortices 
interact with a random distribution of attractive pinning 
centers with U p (r) — —A p e~( r / TpS> ~. Length is normal- 
ized by r p , energy by A v , and time by r = r\r^lA v . We 
consider N v vortices and N p pinning centers in a rect- 
angular box of size L x x L y . Moving vortices induce an 
electric field E = f V x z, with V = £\ dRJdt. 

To study the fluctuation-dissipation relation (FDR), 
we proceed in a similar way as in previous simulations of 
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structural glasses H and consider the observables: 



N v 



N v 



Mt) = TF L i ^(*) = > . . (1) 
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where s 

SiSj 



Sij, and rf 



1, 1 are random numbers with = and 
i?f — i?^ m with fj, = x,y and R cm 
being the center of mass coordinate. Taking F = Fy 
we study separately the FDR in the transverse and par- 
allel directions with respect to F. The time correlation 
function between the observables and is 



t ) = (AJt)B^(t )) 



1 



N v 

W Y,(r?(t)r?(to)) , (2) 

1= 1 



are independent of the Si. The integrated 
response function Xp, f° r the observable is obtained by 
applying a perturbative force f- = e-Sip, (where jj, = x, y) 
at time to and keeping it constant for all subsequent times 
on each vortex: 



1 



X M (t,t ) = lim-[(^ I (t)) e 



(3) 



We then define a function, T e ^ f (t,t ), by the relation, 



2k B T£ t (t,t ) 



A M (t,t ) , 



(4) 



where A p (t, to) = £ Ef^lrf (t)-rf (i )| 2 ) = C M (t,t) + 
C M (to, to) - 2C M (t, to) is the quadratic mean displacement 
in the direction of /i. For a system in equilibrium at 
temperature T the FDT requires that T* it = T e y ff = T. 
In a noncquilibrium system, like the driven vortex lattice 
with pinning, the FDT does not apply. Since we are 
interested in the stationary states reached by the driven 
vortex lattice, where aging effects are stopped pL H, [l8[ , 
then all observables depend only on the difference t — 
to, if we choose to long enough to ensure stationarity. 
From the parametric plot of x M (t) vs. A M (f) one defines 
the effective temperature T£ fi (t) using Eq. (4) , provided 
Tg lt (t) is a constant in each time-scale 

We study the transverse and longitudinal FDR for the 
moving vortex lattice as a function of driving force, F, 
solving the dynamic equations for different values of A p , 
n v , and T. The simulations are performed with pin- 
ning density n p = N p r p / L x L y = 0.14 in a box with 
L x /L y = V3/2 and N v = 256. We consider A p /A v = 
0.35,0.2,0.25,0.1, n v = N v r 2 p /L x L y = 0.05,0.07, and 
T < 0.01. We impose periodic boundary conditions with 
the algorithm of Ref . . Averages are evaluated during 
80000 steps of At = 0.1t after 65536 steps for reaching 
stationarity. To calculate the response function X^t), 
given by Eq. (3), we simulate two replicas of the sys- 
tem, with the perturbative force ff = esip, applied to 
one of them. Starting from the same initial condition 
we let the perturbed and unperturbed system evolve for 
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FIG. 1: Averaged transverse quadratic mean displacement 
A x (t) (a), and transverse response function x^(t) (b), for 
A p = 0.2, T — 0.01 and V = 0.11, in the smectic flow regime, 
(c) Transverse fluctuation dissipation relation Xx(t) vs A x (t) 
for T = 0.01. The dashed line indicates A X /2T, with T = 0.01 
the thermal bath temperature. T e u(T) is obtained from the 
linear fit to x^(t) for A x > (dashed dotted). We find that 
TSf(T) » T£ f (0) + T, where T£ f (0) is the corresponding 
transverse effective temperature for T — (inset). 



5000 time steps and calculate A /J (t) c and A Pi (t) e=Q re- 
spectively. The replicas then evolve again after changing 
the realization of the random factors, Si, and taking the 
final configuration of the unperturbed system as the new 
common initial condition. From this we get (-A^(f)) , 
both for e = and e ^ 0, and thereby the response func- 
tion, x^Wj i s determined. We choose e small enough in 
each case to assure a linear response of (A^(t)) e . 

At T = 0, there are three different dynamical regimes 
|l6f when increasing F above the critical depinning force 
F c : two fluid phases with plastic flow for F c < F < F p 
and smectic flow for F p < F < F t , and a transverse 
solid for F > F t . For fixed pinning density n p , the char- 
acteristic forces F c , F p and F t depend on the disorder 
strength A p and vortex density n v . We start by ana- 
lyzing the FDR in the smectic flow regime for different 
values of A p , n v and T. In Fig. 1(a) we show the typical 
transverse quadratic mean displacements, A x (t), and in 
Fig. 1(b) the integrated transverse response, Xx(t), for 
this dynamical regime. In Fig. 1(c) we show the FDR 
parametric plot of Xx(t) against A x (t). We see that the 
equilibrium FDT does not apply in general but two ap- 
proximate linear relations exist for A x (t) < 0.05?"^ and 
for A x (t) > rz, with a non- linear crossover between them. 
Following Eq. (4), we find that the short displacements 
region corresponds to the bath temperature T = 0.01, 
and therefore the equilibrium FDT applies in the trans- 
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FIG. 2: Averaged longitudinal quadratic mean displacement 
A y (t) (a), and longitudinal response function Xvif) 0>)i f° r 
A p — 0.2, T = 0.01 and V = 0.11, in the smectic flow regime, 
(c) Longitudinal fluctuation dissipation relation R y (t) vs 
A y (t) for T = 0.01. The dashed line indicate A X /2T, with 
T — 0.01 the thermal bath temperature. 



verse direction only for short times, t <C r p /v. For 
the large displacements region we get an effective trans- 
verse temperature T'£ fi (T) = 0.045 > T. fn the inset of 
Fig. 1(a) we show the FDR for T = 0. Comparing the 
results for different T, we find T* lf {T) w T* ti (0) + T. 
In Fig. 2 we analyze the FDR for the longitudinal direc- 
tion: Fig. 2(a) shows the quadratic mean displacement, 
A y , and Fig. 2(b) shows the response, Xy I n Fig- 2(c) 
we obtain the corresponding FDR. We observe that the 
equilibrium FDT applies for A x (t) < 0.05^ at the bath 
temperature T. The plot does not have a constant slope 
for larger displacements. Similar results are obtained in 
the plastic flow regime. 

In Fig. 3 we show the calculated transverse effective 
temperature, T* ff , for T = as a function of voltage (i.e., 
average velocity, V). We observe that above the critical 
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FIG. 3: Transverse effective temperature T^±± vs voltage V 
for A p — 0.35, T — and n v = 0.07, using the diffusion 
relation (filled diamonds) and a generalized Kubo formula 
(open diamonds). Inset: Transverse diffusion constant D x vs 
V . Dashed lines indicate the transverse freezing transition at 
F — Ft and the dash-dotted line indicate the melting temper- 
ature of the unpinned system T m « 0.007. 
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FIG. 4: Transverse effective temperature T^ tt vs voltage V, for 
different values of pinning amplitude A p , and vortex density 
n v , at T = 0. Ap = 0.35, n v = 0.05 (■), A p = 0.35, n v = 0.07 
(□), A p = 0.2, n v = 0.05 (▲), A p = 0.2, n v = 0.07 (a), 
Ap = 0.1, n v = 0.05 (0), A p = 0.1, n„ = 0.07 (♦). The inset 
shows 23* vs V/AI-", with a = 0.5 



force, T% ti , is a decreasing function of V that reaches amplitude, A p , and vortex density, 



a value close to the equilibrium melting temperature of 
the unpinned system, T m 0.007 po|] , when the system 
approaches the transverse freezing transition at F = Ft 
(obtained from the vanishing of the transverse diffusion 
D x , shown in the inset). It becomes very difficult to 
compute numerically T% ti for driving forces F > F t , since 
A x and Xx are bounded at T = while for finite T there 
are very long relaxation times involved. Therefore, we 
leave the interesting case of obtaining T eit (T) for F > 
F t for future study. The opposite limit F — > is also 
interesting. We shall report on FDT measurements when 
F = 0, in the Bragg and vortex glasses, elsewhere. 
In Fig. 4 we show the dependence of T£ f f with pinning 



we observe that T* fi 



T m when F 



In all the cases 
> F t ; even when 



F t depends on A p ,n v . The inset of Fig. 4 shows that we 
can approximately collapse all the curves plotting T® ti 



vs V/Ai 



ith 



0.5. The shaking temperature for 



a single vortex is expected to satisfy this scaling with 
a = |l3|]. The same result should apply in the limit of 
non-interacting vortices or incoherent motion |l . In the 
opposite case of motion of a rigid lattice we can apply 
the one particle result to a Larkin-Ovchinikov correla- 
tion volume where the pinning force summation gives an 
effective pinning amplitude yfA~p, and therefore a = 1 
in this limit. It is noteworthy that the value we find is 
intermediate between these two limits. 
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We now show that the same T* tf can be obtained from 
a different observable. In particular from experimentally 
accessible quantities like the transverse resistivity and 
the voltage fluctuations. If T£ ff is well defined in a given 
time scale we can expect a generalized Kubo formula to 
hold iBLTl, 
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l 



dt'(V x (t)V x {t')) , 



(5) 



where R x (t) = (^^) e=0 - (^^) c=0 is the linear 
transverse resistance. If we now make a parametric plot 
of R x (t) vs. L dt 1 (y x {t)V x {t')) we again find a linear 
slope equal to 1/T for t <C r p /v and a second linear slope 
of l/Tg if for t > r p /v. In Fig. 2(a) we compare the two 
effective transverse temperatures obtained using Eqs. (4) 
and (5). We see that they are similar within the error 
bars. The shaking temperature, T sh , defined in Eq. (3) 
of Ref. H is proportional to the time integral of the cor- 
relation function of the pinning force F pln = fi p (t). 
T sh can be obtained from Eq. (6) if we replace R x (t) with 
the single vortex value Rq = 1/rj and the integral of the 
V x (t) correlation function is taken for all t [since for the 
transverse direction V x (t) cx F* in (t)]. In other words, 
T S h of U corresponds to taking the average slope in the 
parametric plot of the generalized Kubo formula [or in 
the parametric plot of the FDR shown in Fig. 1(c)], see 
also [^TJ. The approach followed here permits defining 
an effective temperature which takes into account all the 
information on its time scale dependence that allows for 
a thermodynamic interpretation of T eff . In this way we 
see clearly that there is a non-trivial value of the trans- 
verse effective temperature T% it for time scales t > r p /v. 
Since one drives the system in the longitudinal direc- 
tion y, the condition of slow motion needed to prove the 
thermodynamic properties in T e y ff is not matched. Thus, 
it is no surprise that we do not find the same value of 
T e ff in the longitudinal direction |], 22 1. Furthermore, 
we find that the short-range correlations of the moving 
fluid are important and give a non-trivial dependence 
with disorder strength A p . Besides this, we have demon- 
strated that in the moving fluid phase T* ff satisfies two 
important results of 0: (i) additivity of temperatures, 
T* it (T) — T* {1 (0) + T and (ii) dynamic freezing occurs 
when T* lf (T) = T m . From the analysis of Eq. (6), we 
note that (V x (t)V x (0)) can be obtained from transverse 
voltage noise measurements (from |V4(w)| 2 ) and R x (t) 
from time dependent measurements of the transverse re- 
sistivity. It will therefore be interesting to have exper- 
imental measurements of T eff to test quantitatively the 
dynamic freezing transition. Finally, we stress that a 
complete dynamic theory of the moving vortex system 
has to capture the features here described. 
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